Entanglement dynamics of one-dimensional driven spin systems 

in time-varying magnetic fields 

Bedoor Alkurtass 1 , Gehad Sadiek 1,2 '*, Sabre Kais 3 

1 Department of Physics, 
King Saud University, 

Riyadh 11451, Saudi Arabia 

2 Department of Physics, 
Ain Shams University, 

Cairo 11566, Egypt 
3 Department of Chemistry and Birck Nanotechnology center, 
Purdue University, West Lafayette, 
Indiana 47907, USA 



* Corresponding author: gehad@ksu.edu. sa 



1 



Abstract 

We study the dynamics of entanglement for a one-dimensional spin chain with a nearest neigh- 
bor time-dependent Heisenberg coupling J(t) between the spins in presence of a time-dependent 
external magnetic field h{t) at zero and finite temperatures. We consider different forms of time 
dependence for the coupling and magnetic field; exponential, hyperbolic and periodic. We exam- 
ined the system size effect on the entanglement asymptotic value. It was found that for a small 
system size the entanglement starts to fluctuate within a short period of time after applying the 
time dependent coupling. The period of time increases as the system size increases and disappears 
completely as the size goes to infinity. We also found that when J(t) is periodic the entanglement 
shows a periodic behavior with the same period, which disappears upon applying periodic magnetic 
field with the same frequency. Solving the particular case where J(t) and h(t) are proportional ex- 
actly, we showed that the asymptotic value of entanglement depends only on the initial conditions 
regardless of the form of J(t) and h(t) applied at t > 0. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 
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I. INTRODUCTION 



Quantum entanglement represents one of the corner stones of the quantum mechanics 
theory with no classical analog Quantum entanglement is a nonlocal correlation be- 
tween two (or more) quantum systems such that the description of their states has to be 
done with reference to each other even if they are spatially well separated. Understand- 
ing and quantifying entanglement may provide an answer for many questions regarding the 
behavior of the many body quantum systems. Particularly, entanglement is considered as 
the physical property responsible for the long-range quantum correlations accompanying a 
quantum phase transition in many-body systems at zero temperature 2|-|4|. Entanglement 
plays a crucial role in many fields of modern physics, particularly, quantum teleportation, 
quantum cryptography and quantum computing 5, (|. It is considered as the physical basis 
for manipulating linear superpositions of quantum states to implement the different pro- 
posed quantum computing algorithms. Different physical systems have been proposed as 
promising candidates for the future quantum computing technology [7 15]. It is a major 
task in each one of these considered systems to find a controllable mechanism to form and 
coherently manipulate the entanglement between a two-qubit system, creating an efficient 
quantum computing gate. The coherent manipulation of entangled states has been observed 



in different systems such as isolated trapped ions [16], superconducting junctions and 
coupled quantum dots where the coupling mechanism in the latter system is the Heisenberg 
exchange interaction between electron spins 18M20I]. One of the most interesting proposals 
for creating a control 



D. Loss et al. 
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able mechanisms in coupled quantum dot systems was introduced by 
221 ] . The coupling mechanism is a time-dependent exchange interaction 
between the two valence spins on a doubled quantum dot system, which can be pulsed over 
definite intervals resulting a swap gate. This control can be achieved by raising and lowering 
the potential barrier between the two dots through controllable gate voltage. In a previous 
work, a two-atom system with time dependent coupling was studied and the critical de- 
pendence of the entanglement and variance squeezing on the strength and frequency of the 
coupling was demonstrated [23 ]. 

Quantifying entanglement in the quantum states of multiparticle systems is in the focus 
of interest in the field of quantum information. However, quantum entanglement is very 
fragile due to the induced decoherence caused by the inevitable coupling to the environ- 
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ment. Decoherence is considered as one of the main obstacles toward realizing an effective 



quantum computing system 



24j. The main effect of decoherence is to randomize the rela- 



tive phases of the possible states of the considered system. Quantum error correction [25 1 
and decoherence free subspace 0, 0] have been proposed to protect the quantum property 
during the computation process. Nevertheless, offering a potentially ideal protection against 
environmentally induced decoherence is a difficult task. Moreover, a spin-pair entanglement 
is a reasonable measure for decoherence between the considered two-spin system and its envi- 
ronment constituted by the rest of spins on the chain. The coupling between the system and 
its environment leads to decoherence in the system and sweeping out entanglement between 
the two spins. Therefore, monitoring the entanglement dynamics in the considered system 
helps us to understand the behavior of the decoherence between the considered two spins 
and their environment. Particularly, the effect of the environment size on the coherence of 
quantum states of the system can be considered by watching the spin pair entanglement 
evolution versus the the number of sites N in the chain. 

Developing new experimental techniques enabled the generation and control of multipar- 



ticle entanglement [28M33I] as well as the fabrication of one dimensional spin chains 
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36]. 



This progress in the experimental arena sparked an intensive theoretical research over the 
multiparticle systems and particularly the one dimensional spin chains [37-46]. The dynam- 
ics of entanglement in an XY and Ising spin chains has been studied considering a constant 
nearest neighbor exchange interaction, in presence of a time varying magnetic field repre- 
sented by a step, exponential and sinusoidal functions of time 47], |48|]. Furthermore, the 
dynamics of entanglement in a one dimensional Ising spin chain at zero temperature was 
investigated numerically where the number of spins was seven at most 49j. The generation 
and transportation of the entanglement through the chain, which irradiated by a weak res- 
onant field under the effect of an external magnetic field were investigated. Recently, the 
entanglement in anisotropic XY model with a small number of spins, with a time dependent 
nearest neighbor coupling at zero temperature was studied too [50]. The time-dependent 
spin-spin coupling was represented by a dc part and a sinusoidal ac part. It was observed 
that there is an entanglement resonance through the chain whenever the ac coupling fre- 
quency is matching the Zeeman splitting. Very recently, we have studied the time evolution 
of entanglement in a one dimensional spin chain in presence of a time dependent magnetic 
field h(t) considering a time dependent coupling parameter J(t) where both h(t) and J(t) 
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were assumed to be of a step function form [5l| . Solving the problem exactly, we found that 



the system undergoes a nonergodic behavior. At zero temperature we found that the asymp- 
totic value of the entanglement depends only on the ratio A = J/h. However, at nonzero 
temperatures it depends on the individual values of h and J. Also we have demonstrated 
that the quantum effects dominate within certain regions of the temperature-A space that 
vary significantly depending on the degree of the anisotropy of the system. 

In this work, we investigate the time evolution of quantum entanglement in a one di- 
mensional XY spin chain system coupled through nearest neighbor interaction under the 
effect of an external magnetic field at zero and finite temperature. We consider both time- 
dependent nearest neighbor Heisenberg coupling J(t) between the spins on the chain and 
magnetic field h(t), where the function forms are exponential, periodic and hyperbolic in 
time. 

This paper is organized as follows. In Sec. II, we present our model and discuss the 
numerical solution for the the XY spin chain for a general form of the coupling and magnetic 
field. Then, we present an exact solution for the system for the special case J(t) = Xh(t), 
where A is a constant. In Sec. Ill, we evaluate the entanglement using the magnetization 
and the spin-spin correlation functions of the system. We present our results and discuss 
them in sec. IV. Finally, in Sec. V we conclude and discuss future directions. 

II. THE TIME DEPENDENT XY MODEL 
A. A Numerical Solution 

In this section, we present a numerical solution for the XY model of a spin chain with N 
sites in the presence of a time- dependent external magnetic field h(t). We consider a time- 
dependent coupling J(t) between the nearest neighbor spins on the chain. The Hamiltonian 
for such a system is given by 

J(t) N Jit) N N 

H = --f-(i+i)E - -7)E - E Kttf , (1) 

z i=\ z i=i i=i 

where <jj's are the Pauli matrices and 7 is the anisotropy parameter. For simplicity, we'll 

consider h = 1 throughout this paper. Defining the raising and lowering operators aj, 

4 = 1(0? + ^), a, = - ■ (2) 
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Follo wing the standard procedure to treat the Hamiltonian (TO) , we introduce Fermi operators 



52) 



,t _ i,t 



i-l 



blexp(i7rJ2b}bj] 



i-l 



a-i = exp(— ITT 



Y.b)b 3 )k 



(3) 



i=i 



then applying Fourier transformation we obtain 



6j 



AT/2 



N/2 



p=-N/2 



J2 ' ''"'''<'»■ 

p=-N/2 



(4) 



where (fi p = Therefore, the Hamiltonian can be written as 

N/2 

H =^2 H p , 

P =i 



(5) 



with H„ given by 



H p = ap(t)[cLcp + cl p c_ p } + iJ(t)S p [clcL p + c p c^ p ] + 2h(t) 



(6) 



where a p {t) = —2J(t) cos0 p — 2h(t) and S p = 2^sm(j) p . 

As [Hi,H m ] = for l,m — 0, 1, 2, . . . , N/2, the Hamiltonian in the 2 Ar -dimensional 
Hilbert space can be decomposed into N/2 non-commuting sub-Hamiltonians, each in a 
4-dimensional independent subspace. Using the basis {|0) , c p c_ p |0) , c p |0) , cL p |0)} we ob- 
tain the matrix representation of H p 



H n 



2h(t) -U(t)8p 

iJ(t)5 p -AJ(t)cos(j) p -2h(t) 

-2J(t)cos0 p 



\ 



(7) 



V -2J(t) cos (j) p J 

Initially the system is assumed to be in a thermal equilibrium state and therefore its 



initial density matrix is given by 



P P (0) 



-PH P (0) 



(8) 



where = 1/kT, k is Boltzmann constant and T is the temperature. 

Since the Hamiltonian is decomposable we can find the density matrix at any time t, 
p p (t), for the pth subspace by solving Liouville equation given by 



ipp(t) = [Hp{t),p p {t)} 



(9) 



which gives 

p p (t) = U p (t)p p (0)Ul(t) . (10) 
where U p (t) is time evolution matrix which can be obtained by solving the equation 

iU p {t) = U p {t)H p {t). (11) 

To study the effect of a time- varying coupling parameter J(t) we consider the following 
forms 



Jexptf) = Ji + (J - Ji)e- Kt 
J CO s{t) = Jo ~ Jo cos (Kt) , 



Jsinit) = J - J sin(Kt) 
J\ — Jo 



Jtanhit) — Jq + 



tanh ( K(t --))+! 



(12) 
(13) 
(14) 
(15) 



Note that Eq. f JTT]) gives two systems of coupled differential equations with variable co- 
efficients. Such systems can only be solved numerically which we adopt in this paper. 



B. An Exact Solution for Proportional J and h 



In this section we present an exact solution of the system using a general time-dependent 
coupling J(t) and a magnetic field with the following form: 



J{t) = A h{t) 

where A is a constant. Using Eqs. ([7]), f lTTj) and f fT6l) we obtain 



1 Mil 


Ul2 1 




' Mil 


"If 


t U21 


U22 J 




I U 2 1 


U 2 2 1 \ 



I -tip 

5 P —4 cos (ftp — | 



and 



i U33 = -2 cos (f) p J (t) u 33 , u 44 = m 33 
Equation ( IT71) can be rewritten as 

i tij = J(t) H'uj . 

for j = 1,2, where 

H' = 



i8„ 



-i5 p 4 cos (f) p — I 









( Ujl 1 


j, Uj = 






^ U i2 J 



J(t) 



(16) 



(17) 



(18) 



(19) 



(20) 



Introducing a unitary rotation matrix 



S 



Using S to diagonalize H' we obtain 



cos 9 er' sin 9 
. e -i<t> s i n Q cos Q 



SH'S 



I Q—\ 



Ai 



A, 



Where the angles <fi and 9 were found to be 

= (n + 1)tt, tan 26* 



2 cos P + 4 



2 ' 



where n = 0, ±1, ±2, . . ., therefore 



sin 20 



5 p 2 + (2cos0 p + f) 



cos 2^ 



2 cos P + 4 



5 p 2 + (2cos0 p +f) 



Finding Ai and A2 we get 



Ai = W <5 2 + (2 cos 4> p + -) - 2 cos P , A 2 = — W 5 2 + (2 cos (j) p + — ) — 2 cos 



Now we define t> 3 - = 5ttj and substitute in eq. (IT^j) we get 



z ^ = (SH'S~ l + ^S" 1 ) ^ . 



Hence 



2 l>,- 



Ax 
A 2 



Vj. 



Solving this equation we obtain 



Vl 



ism 9, *•/>''«"' 1 ' '- 



Finally w is given by 



Un = cos 2 9e-^ Jo J ^ dt ' + sin 2 9e~^ J ^ dt ' ., 
u l2 = -i sin 9 cos 9 je^ 1 fo J ^ dt ' - e** £ J ^ dt ' 



u 2 i = -U12 
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where 



u 2 2 = sin 



2 n -iAi /J /(f )df 



+ COS 



-a 2 / Q * 



M33 = «44 = e 



!i cos </> p J* J(t')dt' 



sin 6* 



cos 6* 



52 + (2cos0p + f) - (2cos0 p + f) 



\ 2 v /5 p 2 +(2cos0 p +f) 



\ 



52 + (2cos0 p + |) + (2cos0 p + f) 
2^ + (2oos^ + f) 



(32) 
(33) 

(34) 
(35) 



III. SPIN CORRELATION FUNCTIONS AND ENTANGLEMENT EVALUA- 
TION 



In this section we evaluate different magnetization and the spin-spin correlation functions 
of the XY model, then we evaluate the entanglement in the system. The magnetization in 
the z-direction is defined as 



A' 



1/JV 



N 



j=i p=i 

where M p = c^c v + cL p c_ p — 1. In terms of the density matrix, it is given by 

_ Tr[Mp(t)] _ 1 ^ Tr[M pPp (t)] 
{ Z) Tr[p(t)} Tr[p p {t)\ ' 

The spin correlation functions are defined by 

ox / qx qx \ ny / oy cy \ cz / cz cz \ 

D l,m — \°l J m/ ' °l,m — \°l °ml ' °Lm ~ \°l °ml ' 



which can be written in terms of the fermionic operators as follows 52] 



(36) 



(37) 



(38) 



Sf,m ~ 4 (BiAi + iBi + i . . . A m _iB m _iA m ) 



cy 



:-i) 



l—m 



(AiBi + iAi + i . . . B m ^iA m -iB m ) 



(39) 
(40) 



Sf, m = ^ (AiBiA m B m ) 



(41) 
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where 



Ai = b\ + b h Bi = b\ - bi 



(42) 



Using Wick Theorem 53[, the expressions (I39I)- (14"TI) can be evaluated as pfaffians of the form 

,771—1 Fl,m 

Pi + 1,Z+1 ' ' ' (+1,771— 1 *<>il+l,m, 



SI 



TPf 



Pm— l,m— 1 Qm—1, 

^771-1,1 




(43) 



/ 



\l—m 



Q.V 



-Pf 



o Pi,i+i Q 1,1+1 

Pi+ii+i 



^l,m — ~^Pf 



\ 



Pi, I Ql,m Pi 

Fl,m Gi >m 
P 





Ql,m— 1 Pl,m 
Fl+l,m-l Gi + i 

Ffn—X,m—1 G m —\ m 
Pm—l,m 



\ 



(44) 



(45) 



where 



^/,m = (-B/An) , -Pi,m = (A-B m ) , Qi m = (AiA m ) , G^ m = (BiB m ) . (46) 

To evaluate the entanglement between two quantum systems in the chain we use the 
concurrence which has been shown to be a measure of entanglement 54j]. The concurrence 
C{t) is defined as 

C{p) = max(0, A a - X b - X c - X d ) , (47) 

where the Aj's are the positive square root of the eigenvalues, in a descending order, of the 
matrix R defined by 

(48) 



(49) 



R = \/VppVpi 

and p is the spin-flipped density matrix given by 

p = {<J y ® (T y )p*{(Ty ® o y ) . 
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Knowing that p is symmetrical and real due to the symmetries of the Hamiltonian and 
particularly the global phase flip symmetry, there will be only 6 non-zero distinguished 
matrix elements of p which takes the form 55J 



(50) 



^ Pi,i p M ^ 

P2,2 P2,3 
P2,3 P3,3 
ypi, 4 p 4 ,4 J 

Hence, the roots of the matrix R come out to be A a = t/Pi,iP4,4 + |pi,4|, = v /p2,2P3,3+ |P2,3|, 



A 



and A, 



VP2,2P3,3 - |P2,3 I 

To find the non-zero matrix elements of p, one can use the formula of the expectation 
value of an operator in terms of density matrix (^G^ = Tr(p G)/ Tr(p) along with the 
magnetization eq. f[3"T|) and the spin correlation functions eq. ([3l?j) - (l4"l~j) which give 



1 



Pi,i 

P2,2 
P3,3 
P4,4 = 



-Mf 



1 171 ' ^Lm 



-Mi 
2 ' 



-M z 
2 " 



S, 



Lm 



1 

4 

1 
4 
1 



1 1 

-M z M z — q z A 

2 m 2 m l,m ' ^ ' 

- l -M z - l -M z m + Sl m + \ 

P2,3 — Si m + Sf m , 
Pl,4 — O l m D l m . 



(51) 

(52) 

(53) 

(54) 
(55) 
(56) 



IV. RESULTS AND DISCUSSION 
A. Constant Magnetic Field 

We start with studying the dynamics of the nearest neighbor concurrence C(i,i + 1) for 
the completely anisotropic system, 7 = 1, when the coupling parameter is J exp as well as 
Jtanh and the magnetic field is a constant using the numerical solution. In Figured] we study 
the dynamics of the concurrence with the parameters Jo = 0.5, Ji = 2,h = 1 and different 
values of the transition constant K = 0.1 and 10. We note that the asymptotic value of the 
concurrence depends on K in addition to the coupling parameter and magnetic field. The 
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larger the transition constant is, the lower is the asymptotic value of the entanglement and 
the more rapid decay is. This result demonstrates the non-ergodic behavior of the system, 
where the asymptotic value of the entanglement is different from the one obtained under 
constant coupling J\. In Fig. [2] we study the effect of the system size N on the dynamics 



(a) 





t(J- 1 ) 



(b) 




(d) 1 



0.15 



t(J- 1 ) 



(c) (d) 

FIG. 1: C(i,i + 1) as a function of t with J = 0.5, Ji = 2,h = 1,N = 1000 at kT = and (a) 
J = J exp , K = 0.1 (b) J = J exp , K = 10 ; (c) J = J tanh , K = 0.1- (d) J = J tanh ,K = 10. 



of the concurrence. We select the parameters Jo = 0.5, J\ — 2, h — 1 and K = 1000. We 
note that for all values of N the concurrence reaches an approximately constant value but 
then starts oscillating after some critical time t c , that increases as iV increases, which means 
that the oscillation will disappear as we approach an infinite one- dimensional system. Such 
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oscillations are caused by the spin-wave packet propagation 



48| . We next study the dynamics 




FIG. 2: C(i,i + 1) as a function of t (units of J" 1 ) with J = J e x P ,Jo = 0.5, Ji = 2,h = 1,K = 1000 at 
kT = and N varies from 100 to 300. 



of the nearest neighbor concurrence when the coupling parameter is J cos with different values 
of K, i.e. different frequencies, which is shown in Fig. |3j We first note that C(i, % + 1) shows 
a periodic behavior with the same period of J if). It has been shown in a previous work 
5l| that for the considered system at zero temperature the concurrence depends only on 
the ratio J/h. When J m h the concurrence has a maximum value. While when J » h or 
J « h the concurrence vanishes. In Fig. [3], one can see that when J = J max , C(i,i + 1) 
decreases because large values of J destroy the entanglement, while C(i,i + 1) reaches a 
maximum value when J = J = 0.5. As J(t) vanishes, C(i,i + 1) decreases because of the 
magnetic field domination. In Fig. H]we study the dynamics of nearest neighbor concurrence 
when J = J S i n . As can be seen, C{i,i + 1) shows a periodic behavior with the same period 
as J(t). We note that we get larger values of C(i,i + 1) compared to the previous case 
J — Jcos- This indicates the importance of an initial concurrence to maintain and yield 
high concurrence as time evolves. Comparing our results with the previous results of time 
dependent magnetic field 



similar to its behavior when h = h S i n , where h S i 



48H . we note that the behavior of C(i,i + 1) when J = J cos is 

h (1 — sin (Kt)), and vice versa. 
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(a) 




(b) 
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' '\ l\ l< l\ i\ /> 

1 I 1 1 1 1 ' | 11 /I 
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1 | ' II 1 | \ 1 1 / 1 
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h(t) 

AM 
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40 50 



t(J~ 1 ) 

(c) 

FIG. 3: Dynamics of nearest neighbor concurrence with 7 = 1 for J cos where J = 0.5, h = 1 at kT = 
and (a) K = 0.1 ; (b) K = 0.5 ; (c) K = 1. 
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i i 1 
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f \ i i 



1 1' 



-urn 
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50 100 150 

t(J~ 1 ) 

(c) 

FIG. 4: Dynamics of nearest neighbor concurrence with 7 = 1 for J sin with J = 0.5, /i = 1 at kT = and 
(a) K = 0.1 ; (b) K = 0.5 ; (c) Jf = 1. 
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B. A Time-Dependent Magnetic Field 



In this section we use the exact solution to study the concurrence for four forms of 
coupling parameter J exp , J tan h, Jcos and J sin when J(t) = Xh(t) where A is a constant. We 
have compared the exact solution results with the numerical ones and they have shown 
coincidence. The dynamics of C(i,i + 1) for h(t) = 1 and J = J exp , Jq = 0.5, J\ = 1 with 



(a) , 




(a) 



(b) 1 

* as-- 
o- 

0.2 
0.15 



-J(t) 
h(t) 





(c) 

FIG. 5: Dynamics of nearest neighbor concurrence with 7 = 1 at kT = 0, Jq — 0.5, J\ = 1, K = 0.1 and 
(a) h(t) = 1 ; (b) h(t) = J(t) = J exp ; (c) h(t) = J(t) = J tanh . 



K = 0.1 is explored in Fig. 5(a) Comparing with Fig. 5(b), which shows the dynamics of 
C(i, i + 1) for h(t) = J(t) = Jexp, Jo = 0.5, J\ — 1 and K = 0.1, as one can see the time- 
dependent magnetic field caused the asymptotic value of C(i,i+ 1) to decrease. A similar 
behavior occurs when h(t) = J{t) = Jtanh, Jo = 0.5, J\ — 1 with K = 0.1 as exploited in 
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Fig. 5(c) Figures [6]^a) and (b) show the dynamics of C(i,i + 1) when h(t) = J{t) = J cos 
and h(t) = J{t) = J sin respectively, where Jo = 0.5 and K — 1. As can be noticed the 
concurrence in this case does not show a periodic behavior as it did when h(t) = 1 in Figs. [3] 
andU 

In Fig. El^a) we study the behavior of the asymptotic value of C(i,i + 1) as a function 
of A at different values of the parameters Jo, Ji and K where J(t) = Xh(t). Interestingly, 
the asymptotic value of C(i,i + 1) depends only on the initial conditions not on the form 
or behavior of J(t) at t > 0. This result demonstrates the sensitivity of the concurrence 
evolution to its initial value. Testing the concurrence at non-zero temperatures demonstrates 
that it maintains the same profile but with reduced value with increasing temperature as can 
be concluded from Fig. [7(b) ■ Also the critical value of A at which the concurrence vanishes 
decreases with increasing temperature as can be observed, which is expected as thermal 
fluctuations destroy the entanglement. 
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\: V: \; wtV 
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t(J- 1 ) 



50 



(a) 



(b) 



FIG. 6: Dynamics of nearest neighbor concurrence with 7 = 1 at kT — with Jo = ^0 = 0.5, K = 1 for (a) 
J cos and h cos ; (b) J sin and h sin • 



Finally, in Fig. [8] we study the partially anisotropic system, 7 = 0.5, and the isotropic 

system 7 = with J(0) = 1. We note that the behavior of C(i,i +JJ in this case is 

I I 

similar to the case of constant coupling parameter studied previously [5l|. We also note 
that the behavior depends only on the initial coupling J(0) and not on the form of J(t) 
where different forms have been tested. 
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FIG. 7: The behavior asymptotic value of C(i,i + 1) as a function of A with 7 = 1 at (a) kT = ; (b) 
kT = 0.5,1. 




FIG. 8: The behavior asymptotic value of C(i,i + 1) as a function of A at kT = with (a) 7 = 0.5 ; (b) 
7 = 0. 

V. CONCLUSIONS AND FUTURE DIRECTIONS 

We have studied the dynamics of entanglement in a one-dimensional XY spin chain 
coupled through a time-dependent nearest neighbor coupling and in the presence of a time- 
dependent magnetic field at zero and finite temperatures. We presented a numerical solution 
for the system for general J(t) and h(t) and an exact solution for proportional J(t) and h(t). 
For an exponentially increasing J(t) we found that the asymptotic value of the concurrence 
depends on the exponent transition constant value, which confirms the non-ergodic behavior 
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of the system. For a periodic J(t) we found that the concurrence shows a periodic behavior 
with the same period as J if). On the other hand for both periodic coupling and magnetic 
field with same period, the concurrence loses its periodic behavior. When J{t) = Xh(t) where 
A is a constant we found that the asymptotic value of the concurrence depends only on the 
initial conditions regardless of the form of the coupling parameter or the magnetic field. In 
future, we would like to study the effect of an impurity spin on the entanglement along the 
driven one-dimensional spin chain. It will be also interesting to study the decoherence of 
a spin pair (quantum gate) as a result of coupling to a driven one-dimensional spin chain 
acting as its environment. 
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